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Abstract 

For stochastic growth models in the Kardar-Parisi-Zhang (KPZ) 
class in 1 + 1 dimensions, fluctuations grow as i 1 / 3 during time t and 
the correlation length at a fixed time scales as i 2 / 3 . In this note 
we discuss the scale of time correlations. For a representant of the 
KPZ class, the polynuclear growth model, we show that the space- 
time is non-trivially fibred, having slow directions with decorrelation 
exponent equal to 1 instead of the usual 2/3. These directions are the 
characteristic curves of the PDE associated to the surface's slope. As 
a consequence, previously proven results for space-like paths will hold 
in the whole space-time except along the slow curves. 

1 Introduction 

The well known KPZ equation for the evolution of a non-linear stochastic 
growth of an interface on a one-dimensional substrate, x h-> h(x, t) (x £ M 
the space and tGl the time), is [28] 

= vAh ( Xi t ) + l\(Vh{x, t)) 2 + 7]{X, t). (1.1) 

v reflects the smoothing mechanism due to the surface tension, A 7^ in 
front of the non-linear term is responsible for lateral spread of the surface 
and irreversibility, and rj is a local noise term. (II. ip corresponds to a situation 
with small slope u = Vh. Generically, one has v(u) instead of the non-linear 
term, with v the speed of growth as a function of the slope u, having the 
property v"(u) ^ 0. 
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Figure 1: Illustration of the PNG height and its dynamics. The bold vertical 
piece is a nucleation. The arrows indicate the movements of the steps. A 
Java animation of the PNG dynamics is available at [16]. 



The long time evolution has a deterministic macroscopic behavior, i.e., 
the limit shape 

h ma (0 = Km (1.2) 

is non-random. To see a non-trivial fluctuation behavior one has to focus on 
an appropriate mesoscopic scale. 

For the KPZ class in 1 + 1 dimensions, the correlation exponent is 2/3 and 
the fluctuation exponent is 1/3 [5,6,22], that is, the rescaled height function 
at time t 

t-(M):= t(iH °'"'"- ) f" K + ° r1 " 1 (1.3) 

will converge to a non-trivial process as t — > oo. 

The rescaling fll.3p makes sense as soon as the macroscopic shape is 
smooth in a neighborhood of £. However, depending on the initial condi- 
tions, fll.ip can produce spikes in the macroscopic shape. If one looks at the 
surface gradient u = Vh, the spikes of h corresponds to shocks in u, and it 
is known that the shock position fluctuates on a scale t 1 ^ 2 . For particular 
models, properties of the shocks have been analyzed, but mostly for station- 
ary growth (see [14, 15] and references therein). In this paper we discuss 
spike-free regions of the surface, where the limit shape /i ma is smooth. 

Limit processes at fixed time 

Recent results obtained by analyzing solvable models in the KPZ class, show 
that the limit processes are not as universal as the scaling exponent, but one 
has to divide into a few subclasses. One of the most studied model is the 
polynuclear growth (PNG) model, which we shortly introduce. At time t, the 
surface is described by an integer- valued height function x \— > h(x, t) E Z, 
x E R, t E R+, with steps of size one. It thus consists of up-steps (_T) 
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Figure 2: Snapshots of the height function for (a) the PNG droplet and (b) 
the flat PNG geometries. 

and down-steps (~L), see Figure [TJ The dynamics has a deterministic and a 
stochastic part: 

(i) up- (down-) steps move to the left (right) with unit speed and disappear 
upon colliding, 

(ii) pairs of up- and down- steps (nucleations) are randomly added on the 
surface with some given intensity. The up- and down-steps of the nu- 
cleations then spread out with unit speed according to (i). 

Analytic results have been obtained in the following cases starting with a flat 
initial substrate h(x, 0) = 0, x G R. 

(a) PNG droplet: the nucleations form a Poisson point process in space-time 
with intensity p(x,t) = 2 for \x\ < t and p(x,t) = otherwise. The limit 
shape is then h ma (£) = 2a/[1 — £ 2 ]+- 

(b) flat PNG: the nucleation intensity is p{x,t) = 2 (or any other constant), 
then the limit shape is simply /i ma (£) = 2. 

These are the two typical cases obtained starting with non-random 
smooth initial condition, and, by universality, one expects the following to 
holds. 

(a) Curved macroscopic limit shape {h'^ a ^ in a neighborhood of £). The 
limit process describing the rescaled surface h resc is the Airy 2 process, A2, 



where k v and Kh are some vertical and horizontal scaling coefficients de- 
pending only on the macroscopic position £. The Airy 2 process has been 



lim h resc (u,t) 



K v A 2 {u/n h ), 



(1.4) 
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discovered in the PNG droplet model by Prahofer and Spohn [32], and later 
shown to occur in several other models [11,21,23,25-27,35]. 

(b) Straight macroscopic limit shape (/i^aCO = in a neighborhood of £). 
The limit process is the Airy! process, A±, 

lim h Tesc (u,t) = K v Ax(u/K h ). (1.5) 

t^oo 

The Airyi process was discovered in the totally asymmetric simple exclusion 
process (TASEP) by Sasamoto [34], see also [8,9]. Recently, we proved (11. 5D 
for the PNG model with flat growth too [10]. All these convergence are in 
the sense of finite-dimensional distributions. 

For a precise a definition of the Airy! and Airy 2 processes, we refer to the 
original papers [9,32,34] and the review [18], in which also known properties 
are summarized, e.g., the one-point distribution given by the Tracy- Widom 
distributions of random matrix [36]. It can also happen that curved part of 
the surface becomes straight (smoothly). Then, the surface will be described 
by the transition process studied in [11]. The last typical case is stationary 
growth, see [3,30,33], to which we come back at the end of the paper. 

Limit processes at different times: along space-like paths 

The above results concern the description of the interface at a given large time 
t. However, they do not say anything about the correlations of the surface 
height at different times. First in [12] for the PNG droplet and more recently 
in [10] for flat PNG, results along space-like paths have been obtained. Let us 
remind what it is meant with space-like paths. A path t = ir(x) in space-time 
is then called space-like if |vr'(x)| < 1 for all x. 

Suppose that we want to analyze the surface of the PNG droplet in a 
neighborhood of the space-time position (£T, T), for some £ G (—1,1) and 
T ^> 1. For simplicity, we choose 7r(£) = 1 and set 

x(u) = £T + uT 2/ \ t{u) = Ttt(x(u)/T). (1.6) 

Borodin and Olshanski prove [12] that for space-like paths 

r h(x(u),t(u)) -t{u)h miL {x{u)/t{u)) 

Inn TFYul = K vA 2 {u{l - £tt (f))/«h) (1-7) 

with k v = (1 — i^ 2 ) 1 / 6 and = (1 — £ 2 ) 5 / 6 . In particular, the correlations 
between h(x(0), t(0)) and h(x(u),t(u)) with 7r'(£) = are identical to the 
correlations between h(x(0), t(0)) and h(x(u n ),t(u n )) if u n (l — £tt'(0) = u. 
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Geometrically, these points form a line with constant slope in the (x, t) dia- 
gram (these are the bold lines in Figure [3] (right)). 

For the flat PNG, the parametrization is slightly simpler because of 
translation-invar iance: 

x(u)=uT 2/3 , t(u) = Tn(x(u) /T) (1.8) 

with 7r(0) = 1. In [10] we prove that for space-like paths it, 

^ h^t^ = 2 ys Mu/2V% (LQ) 

Slow space-time curves with decorrelation exponent r — 1 

At first sight surprisingly, the scaling coefficients in (II. 9p do not depend on the 
slope of the space-like path. Similarly, in (11.71) the correlations are the same if 
we follow lines with constant angle. Why? The reason is that the space-time 
is fibred in a non-trivial way: along certain directions the decorrelation is 
much more slow, on a T T scale for some scaling exponent r > 2/3. 

Main result. By considering the PNG model we show that r = 1. 

For the flat PNG, the particular curves are parallel to the time axis, while 
for the PNG droplet are the one keeping a constant angle. In general, these 
curves are such that the gradient of the surface is constant: the characteristics 
curves of the following PDE. Let u = V/i be the slope of the surface. Then, 
on a macroscopic scale, dh/dt = t> eff (u), where v eS is the effective speed 
including the noise effect in (11.11) . Then, taking the spatial derivative, we get 

<9u , ,9u . . <9t> eff (u) 

^ + a(u)- = 0, .(») = — (1-10) 
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The characteristics of this equation are the trajectories satisfying dx/dt = 
a(u) and du/dt = 0. 

For the polynuclear growth model, one has + u 2 [30] , from 

which one verifies that the slow directions of Figure [3] are indeed the char- 
acteristics. Why are precisely the characteristics curves to have slow decor- 
relations? Because the information flows along them, so on a macroscopic 
scale, the fluctuations of the height is influenced by the randomness occurring 
around the characteristics. 

In general, since characteristics can cross [13] and one has the phenomena 
of shock formation. In that case, the derivative of the limit shape has a 
discontinuity. As said in the introduction, in this paper we focus away from 
this situation. In terms of the slope u, we consider cases without shocks, i.e., 
either with constant u or in the rarefaction regions. 

We believe that the same phenomena occurs in higher dimensions too. In 
a recent work with A. Borodin [7] on a model in the anisotropic KPZ class 
in 2 + 1 dimensions a similar phenomena of slow curves seems to occur. 

Consequences 

A first consequence is that, as soon as we do not follow the characteristics, 
( 11.91) and (jl.7p hold without the space-like constraint. To illustrate it, let us 
state the precise result for joint distributions in the PNG droplet, which is 
proven in Section 12.31 

Theorem 1. Consider the parametrization $1.6\) but this time tt is any space- 
time path with 7r(£) = 1 and 7r'(£) ^ l/£, i.e., we exclude the case of the 
path tt following the characteristics. Then, in the sense of finite- dimensional 
distributions, 

r h(x(u),t(u)) -t(u)h ma (x(u)/t(u)) 

km — -YT^ = k v A 2 {u(1 - £vr £ ) K h ) 1.11 

with k v = (1- ^ 2 ) 1 /6 an d Kh = (1 - i 2 f/ & . 

The analogue for the flat PNG is that fll.9p holds for any finite tt', since 
in that case the characteristics corresponds to tt' = oo. Therefore, to get a 
different behavior from the Airy processes for the height-height correlation, 
one has to rescale the component of the space-time path along characteristics 
with the scaling exponent r = 1 and the other component using the scaling 
exponent 2/3. For example, for the flat PNG, one has consider the scaling 

x(u)=uT 2/3 , t(v)=T(l + v). (1.12) 
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Then, the observable to be analyzed is the two-parameter limit process of 
the height function rescaled as 



,resc/ S _ h{x{u),t(v)) -t(v)h mSu {x{u)/t(v)) 

The description of this limit process remains an open problem. 



Outline 

The rest of the paper is organized as follows. First we give a heuristic argu- 
ment for r can not be strictly less than one. Then we prove it for the PNG 
droplet (Theorem [2]). This needs also a technical control of the maximum 
(Corollary [U]) for which tightness is needed. For the flat PNG geometry a 
similar result will hold (Claim [7]). In these two cases, one actually sees that 
the heuristic argument is the right interpretation. Finally, in Proposition [8J 
we show the same result for stationary growth, which turns out to be much 
simpler, but in the proof there is no trace of the mechanism. 
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2 Results and proofs 

In this section we prove that for the PNG model, for any r < 1, if we take 
two points along the special directions at time T and time T + T r , then they 
are trivially correlated on the fluctuation scale T 1 / 3 . But first let us present 
a heuristic argument, which then it is transformed into a proof for the PNG 
model. 

2.1 Heuristic argument 

On a macroscopic level, the information flows along the characteristic curves. 
Consider two positions on one of such curves, say A at time T and B at time 
T + T T (the white and black dots in Figure H]). We have the following facts: 

(1) The correlation exponent is 2/3, therefore the fluctuations of the 
heights at A and B are correlated with the randomness occurring in a 
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Figure 4: The white and black dots are two space-time positions lying on 
the same characteristic curve (the bold curve). The height fluctuations along 
this curve are correlated with events in space-time inside a small strip of 
order T 2 / 3 . In the zoom, it is indicated that the correlations of the events 
after time T lies in a smaller region, of width at most T" , with v < 2/3 for 
t < 1. 



strip of order T 2 / 3 around the characteristics. However, if r < 1, the 
events occurring in the time interval [T, T + T T ] will be relevant only 
on a smaller region, of order T v with v = 2r/3. 

(2) Again, because the correlation exponent is 2/3, the fluctuations of the 
height function at time T in a ^-neighborhood of A differ from the 
fluctuations of the height at A by 0(T U ' 2 ) = o(T 1 / 3 ). 

(3) The fluctuation exponent is 1/3, meaning that during a time interval 
T, fluctuations are created on a T 1 / 3 scale. 

By (1), the height at B is the sum of the height at some point in a T v - 
neighborhood of A plus the contribution during the remaining T T time in- 
terval. By (2), the fluctuations of the first term is equal to the fluctuations 
of the height at A up to an o(T 1 / 3 ) term. By (3), the fluctuations in the 
second term are only of order T r / 3 = o(T 1 / 3 ). Therefore, on the natural 
scale of fluctuations, T 1//3 , the height of A and B are trivially correlated in 
the T — > oo limit (i.e., the difference is only a deterministic factor). 

Remark that, like in the PNG model defined below, often the characteris- 
tics curves are straight lines. However, for generic time-dependent dynamics, 
this is not the case. 
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Figure 5: Height and directed polymers for the droplet geometry. 



2.2 PNG and directed polymers 

In the proof, we use the picture of directed percolation on Poisson points, 
whose connection to the PNG model is the following. Let us start with the 
PNG droplet case. The nucleations in the PNG model form a Poisson point 
process in space-time with intensity p(x,t) = 21(|x| < t). We turn by 7r/4 
the point of view as in Figure El To avoid annoying v2 factors, we then 
zoom out by a linear factor \/2, so that the density of Poisson points in the 
directed percolation picture is 1. 

Consider all the directed paths consisting of a sequences of Poisson points, 
where the points are joined by segments of slope in [0, oo]. The length of a 
path is the number of Poisson points visited. Then, for a point (u, v) G 
let L((0, 0), (u, v)) be the length of the longest path from (0, 0) to (u, v). We 
call maximizers the paths of longest length. As noticed in [31], the height 
function of the PNG droplet at position (x,t), h(x,t), can be expressed by 
L, namely 

h(x,t) = L((0,0),(x + t,t-x)). (2.1) 

For the flat PNG, the only difference is p(x, t) — 21(t > 0), which implies 
that the height function is given by (12.11) but for L the optimization is made 
line-to-point: from {u + v = 0} to (x + t,t — x). More details can be found 
in the literature, for example in Section 3 of [19]. 

2.3 PNG droplet 

Denote by h(x,t) = h(x,t) — th ma (x/t) the difference between the height 
function and the approximate position from the limit shape. 
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Theorem 2 (PNG droplet). Consider a fixed fixed £ G (—1,1). For any 

r < 1, there exists a [3 < 1/3 such that 

lim P(|h(£T,T) - h(£(T + T T ),T + T T )\ < T 13 ) = 1. (2.2) 

T — >oo 

To transform the heuristic argument into the proof of Theorem [2} we use 
results of [21] which control the geometric distribution of the maximizers 
close to the end points. The results in [21] use ideas from in [24] about 
transversal fluctuations. Then, it will only remain one technical difficulty: 
the convergence to the Airy2 process of [12, 32] are in the sense of finite- 
dimensional distributions. This is not enough to conclude that the maximum 
converges: a tightness result is needed. 

Proof of Theorem^ First, notice that it is enough to consider the case £ = 0. 
Indeed, if we consider a point (£T, T) and apply the mapping 

$ . ( X; t ) _> (i _ £ 2 )-V2 (x -t£,t- ^), (2.3) 

we recover the situation with £ = with time given by a/1 — £ 2 T instead 
of T. This mapping preserves the Poisson measure on nucleations and, in 
the percolation picture directed paths remain directed. Therefore, w.l.o.g, 
we consider the £ = situation. 

As illustrated in Figure EJ denote O = (0,0), A = (0,T), and 
B — (0, T + T r ). Moreover, let C(r, v) be the cylinder of axis AB and width 
w = 2T U , and dC(r, v) be the border border of the cylinder without lids. We 
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consider the following event 

D = {Poisson points configurations such that the maximizers 

from O to B do not cross dC(r, u)}. (2.4) 

This quantity was studied in the proof of Theorem l(ii) of [21]. To be more 
precise, one uses the bound on P(f2 \ D) in the equation just after (20) with 
the following improvements. The e is the one of Lemma 2 and it is actually 
t-dependent by using (11). The bound on P(-F) is in the equation after (22), 
which is also t-dependent by using the bound of [2] reported in (8). In the 
case r < 1 and setting v — | — we obtain 

P(D) > 1 - ci exp (-c 2 T^- T ^ 12 ) (2.5) 

for some constant C\, c 2 > independent of t. 

Consider the hyperbole 7i, along which the height is macroscopically con- 
stant (~ 2T): 

H = {(x,t) = (r(6) sm(6),r(9) cos(9)),r(9) = Ty/cos(20), 6 e (-7r/4, tt/4)} 

(2.6) 

and let S(u) be the piece of Ti lying inside the cylinder C(r, v). 6 in (12.61) is 
the angle between the x = axis and the segment from O to the point in Ti. 

The bound (I2.5P tell us that, with probability going to one as T — » oo, 
the maximizers of L(0,B) crosses 5'(z/). Therefore we need to concentrate 
on those situations only, for which we have 

L(0, A) + L(A, B) < L(0, B) < L(0, S) + L(S, B). (2.7) 

We rewrite it as 

L(A,B)-2T T < L(0,B)-L(0,A)-2T T < L(0, S)-L(0, A)+L(S, B)-2T T . 

(2.8) 

Recall that L(0, A) ~ 2T, L(0,B) ~ 2(T + T r ), and L(A, B) ~ 2T T . Thus, 
we have to show that both right and left hand side are (in absolute value) 
smaller than T 13 with probability 1 as T — > oo, for some (5 < 1/3. 
From Lemma 7.1 of [2], setting (5 — (1 + r)/6, we have 

P{\L{A, B) -2T T \> T 13 ) < ci exp(-c 2 T (1 - T)/4 ). (2.9) 

Moreover, since L(A,B) < L(S,B) < L((x,t = T),B) and both terms are 
2T T + 0(T T / 3 ) [3,31], by Markov inequality it follows 

P(\L(S, B) — 2T T \ > T 13 ) < 0(1/T^ T ^ 6 ) —>■ as T ^ oo. (2.10) 
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A better bound of the type (12.91) can be also obtained using the results of [4], 
but it is not needed for our purposes. 
Finally, consider the points 

X{u) = (x{u),t(u)) on H (2.11) 

associated with the angle 9 = uT~ 1/3 . Then, u i-> (L(0,X(u)) - 2T)/T 1/3 
converges, in the sense of finite-dimensional distributions, to the Airy2 pro- 
cess Aziu) [12,20]. The Airy 2 process is a.s. continuous and the points in 
S(u) corresponds to \u\ < T u ~ 2 ^ 3 = 1/T^ 1_T ^ 6 . The control of the maxi- 
mum requires a bit more than finite-dimensional distributions, which is the 
content of the next part of the paper leading to Corollary [6j there exists a 
P < 1/3 such that P(\L(0,S) - L(Q,A)\ > T p ) -> as T -> oo. This end 
the proof. □ 

With this Theorem we can then prove Theorem [TJ 

Proof of Theorem^ Using the transformation (12. 3p . we consider w.l.o.g. the 
£ = case. As in Theorem [21 denote h(x,t) = h(x,t) — th ma (x/t) and 
consider m points parametrized by space-time positions Xi = UiT 2 ^ 3 and 
ti = T-K(uiT~ 1 / 3 ). Then we want to determine the large T limit of 



P 



f){h{x k ,t k ) < s k T 1 ' 3 } = P f]{h(x hl T) < s k T 1 ' 3 + E k } ) (2.12) 
\fc=i / \fe=i 



where S k = h(x k ,T) — h(x k ,t k ). Denote the events Q k = {|Hfc| < T^}, and 
A k = {h(x k , T) < SfcT 1 / 3 + E k }. Then we can rewrite (12.121) as 

(EH) = p (f) {{A k n Q k ) u (A k n QD) \ = P ( (f]( A k n Q*) j u i? 

(2.13) 

where in i? are terms with at least one Aj fl Q c y Then, using inclusion- 
exclusion principle, we get 

tm \ / m \ 

f|(A fc n Q fc ) J + P(i?) - P I R n fj(A fc n Q fc ) j . (2.14) 

By Theorem [21 there exists a. (3 < 1/3 such that P(Qj) — >• 1 as T — > oo, i.e., 
P(Qp — *■ ^" Since the last two terms contain at least one they go to zero 
as T — > oo. For the first term in the r.h.s. of (12.141) . because (3 < 1/3, we 
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Figure 7: Illustration of the multilayer version of the polynuclear growth 
model in the PNG droplet geometry at time t = T. 



have 



lim P f]((h(x k , T) < s k T l ' z + S fc ) n {|S*| < T p }) 

\fc=i / 

(m \ / m 

f](h(x k ,T) < SkT 1 / 3 )) = P ( f)(A 2 {u k ) < s k )) , (2.15) 

where the last equality is the fixed time result of [32]. □ 

Multilayer PNG droplet 

To obtain the result of Corollary [6j we use the multilayer extension of the 
polynuclear growth model introduced in [32]. Let us just recall it briefly. 
Instead of a single height function h(x,t), we have a set of height functions 
{he(x, t),£ < 0}, where h (x,t) coincides with the original h(x,t). Initially, 
he(x,0) = £, for all x G R. A nucleation at level £, for i < —1, occurs 
whenever there is a merging at level I + 1 and the nucleation is at the same 
position of the merging. Moreover, the deterministic part of the dynamics 
for the PNG applies to the new height functions as well (see (i) at page [3]). 
By construction the lines do not intersect (see Figure [7] for an illustration). 
One then defines an extended points process r\ by 



1, he(x,t) = j for some 
0, otherwise 



< 0, 



(2.16) 



13 



Along the hyperbole TC defined in ( I2.6P the point process r\ is stationary, if 
we choose as parametrization the angle 9 with the origin [12]. For any point 
(xi, ti) G TC, we denote by 9i the corresponding angle. As shown in [12,20], r\ is 
an extended determinantal point process, i.e., it has determinatal correlation 
functions: p (n) (6>i, ji, . . . ; 9 m , j m ) = det[K T (9 k , jk', &i, ji)]i<k,i<m with K T the 
extended Bessel kernel. Define the operator H on £ 2 (Z) by 

(ffV)O') = -V>0' + 1) - 1>U - 1) + 2 V>U) + ~V>0') (2.i7) 

and the one-time Bessel kernel by 

with J n (2T) = jL = i dze T ^ _1 ^^ _n_1 the standard Bessel functions [1]. 
Then, the stationary extended Bessel kernel is given by 

K T (9 1 ,j 1 ;9 2 ,j 2 ) = I ~ (2.19) 

Reseated process. Recall from (12.111) that for the rescaled process on TC the 
angle is 9 = u/T 1 ^ 3 , which is thus given by 

^ sc (u u Si ) = T^rjTimT- 1 / 3 , [2T + ^T 1 / 3 ]). (2.20) 

Then, denote the rescaled height function by 

h resc^ = h(x(u),t(u)) - Th m& (x{u) / t{u)) ^ ^ 

with (x(u),t(u)) the coordinate on TC with angle 9 = ujT 1 !' 3 . Also, for a 
given function /, we define 

f T (x)=f((x-2T)/T 1 / 3 ), (2.22) 

so that 

VT>, f) = ^73 E f( S )VT SC (u, S)=J2 fTti)VT(uT-V 3 ,j). 

sG(Z-2T)/T 1 /3 jeZ 

(2.23) 

After this setting, we can go to the proof. The strategy is almost identical 
as in Section 5 of [26]. The only difference, is the derivation of the first lemma 
(Lemma 5.1 in [26]), which is carried out as in Appendix A of [32], exploiting 
the stationarity and the semi-group structure in the kernel. 
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Lemma 3. Assume that f is a C°° function and there exists a constant M 
such that f(x) = if x < —M and f(x) — 1 if x > M. Then, there is a 
constant C(f) such that 

E {(vT C (ui, f) - Vt s >j, /)) 4 ) < C(f) flu, - Uj \ 2 + \ Ui - u 3 \/T 2 ' 3 ) (2.24) 

for — 1 < Ui,Uj < 1 and for allT > T , with T a fixed large enough constant. 

This lemma is the corresponding of Lemma 5.1 in [26] if one takes into 
account that in [26] time is discrete, which means that \ui — Uj\ > 1/T 2//3 . 

Proof of Lemma\^ First of all, using (I2.23P we get for the l.h.s. of (12.241) . 

e((^ sc (^,/)-^ csc K-,/)) 4 ) 

= E ^ E Mil) ■ ■ ■ /t(u)E (vAu^T^, iO--- vAut^ 3 , k)) 

n=0 ^ ' u,...,i 4 eZ 

(2.25) 

where = itj if k > n and = Uj if k < n. The last term in the 
expected value is just the 4-point correlation functions, which is given by a 
4x4 determinant. 

Now, the reason of considering the height function on the line 7i, is that 
we have a stationary determinantal point process with kernel K?- is 
the evolution of one-time kernel by the operator H . This is exactly the 
same setting of Appendix A of [32]. One can copy step by step the proof of 
Lemma A.l in [32] except the last one, which is the only one which uses the 
particular structure of the operator H. The result, is a bound C(f)\ui — Uj\ 2 
plus the following term (see (A. 16) of [32]) 

(m - Uj )T 2/3 Tr (B T F T ) (2.26) 

where Ft is given by 

F T := [H, f T ]f T - f*[H, f T ] - 3f T [H } f T ]f 2 + 3f 2 [H, f T ]f T . (2.27) 

It is not difficult to get an expression for Ft (just use a matrix representation 
for H and fr)- We obtain 

*<U) ' ! - h I=l< (2.28, 

Therefore, 

Tr (B T F T ) = ]T B T (i, j)F T (j, i) = 2 ]T B T (i, i + l)(f T (i + 1) - / T (z)) 4 , 

(2.29) 
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where we also used the symmetry of the kernel B^, given in (12.181) . Recall 
that we have f T (i + 1) - f T (i) = if i <£ [2T - MT 1 / 3 , 2T + MT 1 ' 3 } n Z. 
Also, \f T (i + 1) - f T (i)\ ~ II/'IU/T 1 / 3 . Thus, 

IHI< const ll/'I^T- 4 / 3 £ |5 T (i,i + l)|. (2.30) 

ie[2T-MT 1 / 3 ,2T+MT 1 / 3 ]nZ 

It remains to see that the last sum is bounded by a constant. The sum is 
bounded by 

^ ^ I^[2T+(s+a)t 1 /3](2T)|| J [2T+(s+A)T i/3 +1 ](2T)|. (2.31) 

By using the uniform bounds (for T > Tq, Tq fixed) reported in Lemma A.l 
of [17] (the first bound is due to Landau [29]): 

(T 1 / 3 J[ 2T+sT i/3](2T)| < const , for all s, 
\T 1/3 J [2T+sT i /3] (2T)\ < const e~ s/2 , for s > 0, (2.32) 

it follows that | (12.311) | < const with the constant independent of T. Putting 
all together, we get 

| fl2T26|) | < const \ Ui - Uj\T~ 2/3 , (2.33) 

with const depending only on /. The uniformity of C(f) for T > T Q in the 
term C(f)\ui — Uj\ 2 is also verified using the exponential bounds (I2.32p . □ 

For a short while, consider only the point process at a discrete set of 
angles 0< = w^T" 1 / 3 , with m = i/T 2 / 3 , ieZD [-T 2 / 3 , T 2 / 3 ]. 

Theorem 4. Let h r ^ sc (u) be defined on u G (T- 2 / 3 Z) n [-1, 1] and by linear 
interpolation. Then, there is a continuous version of A2{u) and 

h^ sc {u) -> A 2 (u) (2.34) 

as T — > oo in the weak* -topology of probability measures on C([— 1, 1]). 

Proof of Theorem^ For Ui = i/T 2 ^ 3 and integer i, the bound of Lemma [S] is 
simply given by C(f)\u{ — uj\ 2 because \ui — Uj\ > T~ 2 ^ 3 . This is the analogue 
of Lemma 5.1 of [26]. Moreover, all the work of Section 5.2 of [26] is valid 
here too, without modifications. Since the convergence of finite dimensional 
distribution was proven in [12], we have weak*-convergence of h™ sc (it is the 
corresponding of Theorem 1.2 of [26]). □ 

The final step is to use Theorem H] together with a control of the maxi- 
mum in the intervals between points in Z/T 2 / 3 to get the convergence of the 
maximum. 
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o x 



Figure 8: Construction for the small distance bound with two points at 
distance 0(1) on H: X\ = X( Ui ) and X 2 = X(m + T~ 2 / 3 ). 



Lemma 5. For i E J = Z D [— T 2 / 3 ,T 2 / 3 ] we have the uniform bound (for 
T >1) 

P ( max h™ sc ( Ul + x)- meQL{h^ c (ui), h^ sc ( Ui + T" 2 / 3 )} > T- 1/6 ] 

\0<x<l/T 2 / 3 ' ' J 

< const e~ Tl/ \ (2.35) 

from which it follows 

P ^ m&x^ sc ( Ui + x) - maxh^iui) > T~ 1 '^ < const T 2 / 3 e' Tl/ \ 

(2.36) 

Proof of Lemma Let 71 be the shaded region of Figure [HJ Then, 

tresc/ \ r;resc/ 1 zresc/ , rr--2/3\l ^ # PoisSOn points in TZ 

max h T rS sc (ui+x)-max{h^ sc (ui),h^ sc (ui+T z/i )} < ^ . 

0<x<l/T 2 / 3 ' T L ' d 

(2.37) 

where the factor T -1 / 3 comes from the rescaling in h r r£ sc . The area of 1Z is 
smaller than 1. Therefore 



l.h.s. of (12351) < P(# Poisson points in ft > T 1/6 ) 

const 
n! 



Econsi T i/6 
— < const e~ T , 2.38 

n>TV6 

for T > 1. □ 

A direct consequence of Theorem H] and Lemma [5] is resumed in the fol- 
lowing corollary. 

Corollary 6. For any 5 > 0, there exists a /3 < 1/3 snc/j t/iat 

lim P f max T 1/3 (h™ sc (u) - ^ osc (0)) > J = 0. (2.39) 
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B 




Figure 9: Geometric construction for the flat PNG. The coordinates of the 
points are O = (0, 0), A = (0, T), and B = (0, T + T T ). The maximization is 
made from the (R, 0)-axis to the points A and B. 



2.4 Flat PNG 

Now that we looked in detail the case of the PNG droplet, we just state the 
corresponding claim for the flat PNG. We present the argument to prove 
it, but we will not do the tightness argument for that case. Recall that 
h(x, t) = h(x, t) - 2t for flat PNG. 

Claim 7 (Flat PNG). For any r < 1, there exists a (5 < 1/3 such that 

lim P(|h(0, T) - h(0,T + T T )\ < T 13 ) = 1. (2.40) 

The strategy is almost identical to the proof of Theorem [5] for the PNG 
droplet. The difference are the following. The line Ti along which one has a 
stationary process is just the straight line (R, T). From the same analysis of 
the transversal fluctuations in [24], we get that for any e > 0, the maximizers 
of I/((R, 0), B) start in a T 2/,3+e -neighborhood of the origin O with probability 
one in the T —>■ oo limit. We can choose for example e = (1 — t)/2, which 
implies that the straight line from (T 2 / 3+e ,0) to B crosses the lines (R, T) 
at a position of order with /x = | — Then, as we made for (12.51) . 
from [21] it follows that ( T — > oo) the maximizers in the time interval 

[T, T + T T ] lies inside a cylinder of width T u , with v = | — The rest of 
the argument is the same, but to transform into a mathematical proof one 
would still need a tightness result, which we do not present in this paper. 
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2.5 Stationary PNG 



Finally we consider the stationary growth of the PNG model. The argument 
is the simplest one, but not very instructive of what actually happens. 

The stationary PNG with mean slope zero is described in Chapter 3 
of [30]. The up- and down-steps of the height function are two independent 
Poisson point processes with density one as shown in Proposition 3.1 of [30]. 
Denote by h stat (x,t) the height function at position x and time t, where the 
reference frame is set at time t = by /i stat (0, 0) = 0. 

Proposition 8 (Stationary PNG). For any r < 1, there exists a (3 < 1/3 
such that 

lim P(|/i stat (0, T + T T ) — /i stat (0, T) -2T T \< T p ) = 1. (2.41) 

T — >oo 

Proof of Proposition 0. By stationarity, 

P(|/i stat (0, T + T r ) -/i stat (0,T) -2T T \ < T' 3 ) = P(|/i stat (0, T T ) — 2T T \ < T 13 ) 

(2.42) 

Set (3 G (t/3, 1/3) and use the notation L = T T and $ = (3/r > 1/3. Then, 
([232D = P(|/^ stat (0, L) - 2L\ < L p ). (2.43) 
From [31] we know that 

lim P(/i stat (0, L)-2L< sL 1/3 ) = F {s), (2.44) 

L^oo 

for some well-defined distribution function F (for details on F see [3,31]). 
Therefore limr^.^ fl2.43p = 1, because (3 > 1/3. □ 
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